One-dimensional quantum wires: 
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1 Introduction 

In these lecture notes we will consider systems in which the motion of electrons is confined 
to one dimension (ID). In these so-called quantum wires electron-electron interaction effects 
play an important role because the restricted dimensions enhance the scattering between the 
electrons and completely destroy the quasi-particle picture. New density wave excitations 
appear that are described by bosonic operators. Here we will develop this bosonic description, 
following a "pedestrian" approach which does not require any previous knowledge in field 
theory methods. These notes therefore serve as a detailed introduction into bosonization by 
carefully deriving the most fundamental formulas. For advanced topics we recommend to 
consult one of the more sophisticated reviews on bosonization [1-5] for further reading. 

1.1 What are quantum wires? 

Before we address many-body effects, let us review some introductory material in order to 
define what quantum wires are, what can be measured in typical experimental setups and 
what theoretical models we wish consider. 



1.1.1 Band structure 

In the classical world we have some intuition on how to create a one-dimensional transport 
channel for particles. In order to make it truly one-dimensional we could simply narrow a pipe 
until two particles can no longer pass each other, for example by creating a narrow path for 
pedestrians until they can only walk in single file. It is easy to imagine how a stop-and-go 
density wave will be created in such a situation. In quantum physics, particles are instead 
described by wave- functions, which are known to be given by discrete standing waves when 
we confine an electron with box-like boundary conditions as for example shown in Fig. Q]for 
a simplified square wire. The allowed standing waves must have wave-numbers that fit in the 
square potential along the x— and y— directions, i.e. d = n x ^f = n xjr- and d = n y ^- = ny^- 
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Figure 1: Confinement in a square wire 

while the motion in the z-direction is unrestricted. The allowed k— vectors are therefore along 
quantization lines which cut the band-structure of the material of the wire. This situation 
is depicted in Fig. [21 where only the x— and z— direction is shown. The original Fermi 
surface of the material of the wire is schematically drawn dotted in this plane. Small energy 
excitations are now only possible along the quantization lines and close to the Fermi surface. 
Each quantization line that crosses the Fermi surface therefore corresponds to a "channel", 
which contributes to the conductivity with quantized conductance as we will see later. 




Figure 2: Cross section of Fermi surface. Vertical lines are the ID quantization lines. Right: 
Effective one-dimensional band-structure along the z-direction. 

For a truly one-dimensional wire we would like to have only one single channel that crosses 
the Fermi surface, i.e. only two Fermi points corresponding to left- and right-moving electrons. 
For this to be true, the energy spacing between the quantization lines must be larger than 
the depth of the Fermi-surface as shown in Fig. [3J For an order-of-magnitude estimate of 
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Figure 3: Electron bands where the energy spacing is larger than the depth of Fermi sea. 



the required length scale of d we can assume that the bandwidth is about a few eV and the 
dispersion is given by an effective mass E = h?k 2 /2m* . Together with k = nir/d the condition 
AE > AE F gives AE = ^ (^f > leV . Using IRy = « 13.6eV and a « 0.5A , we 
finally arrive at the estimate: 



We therefore find that diameters of one nanometer or less are required in order to ob- 
serve one-dimensional physics. For semi-conductors the bandwidth and the effective mass can 
be much smaller, so that a few nanometers may be possible. For metallic wires, however, 
essentially single atomic chains are required. 

1.1.2 Creating quantum wires: Experiments 

The creation of one dimensional wires is therefore still a great challenge even with today's 
semiconductor technology, which can create structures that are only a few tens of nanometers 
in size. However, even if traditional processing techniques of lithography and etching can 
be scaled down to about lOnm, the resulting structures are still imprecise on a nanometer 
scale. The resulting wire typically has relatively uncontrolled wavy edges, which immediately 
leads to localization in a one-dimensional system. In order to produce a wire that has perfect 
translational invariance, we therefore require production mechanisms that give full control 
down to the atomic scale. We will briefly describe three of the most promising experimental 



d < TT 
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Figure 4: ID wire created by cleaved edge overgrowth. [6] 




Figure 5: A carbon nanotube. 



approaches below. 

Using molecular beam epitaxy it is possible to grow atomically perfect structures layer 
by layer. In this way it is now standard to create high quality two-dimensional electron 
gas (2DEG) systems by use of an inversion layer. In the so-called cleaved edge overgrowth 
technique such a sample is then cleaved in-situ [6]. The cleaved edge is typically along one 
of the crystal axis and also atomically perfect. If a gate is overgrown on the edge, it is then 
possible to apply a suitable potential that forces the mobile electrons to be trapped in the 
one dimensional wire that is defined by the intersection of the 2DEG and the cleaved edge as 
shown in Fig. [H A suppression of the density of states which is consistent with the theory has 
been demonstrated in such a setup [7]. 

Another approach to create atomically perfect structures is to use macromolecules that 
form a tubular structure. The most famous examples are the carbon nanotubes as shown 
in Fig. Even though the detailed diameter and location of the tubes cannot be fully 
controlled in the production process, the chemical composition guarantees that the resulting 
tubes are atomically perfect structures on the nanoscale. Several experiments showed that 
single wall carbon nanotubes exhibit the signatures of a one-dimensional many-body system 
[8-10] consistent with theory [11-13]. 
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Figure 6: Atomic Au chains deposited on Si step edges. [14] 



Superstructures on extremely clean surfaces, like a stepped silicon surface are a third 
example [14]. Cleaving a Si(lll) surface at a slight angle creates a regular step structure 
on which small amounts of gold can be deposited. The gold atoms tend to form atomic 
chains along the steps of the silicon, a so-called Si(lll) 5x1 Au structure as shown in Fig. [6J 
Experiments were able to find the signature of spin-charge separation using photoemission on 
such samples [14]. 

Other possible experimental realizations of one-dimensional wires include quantum hall 
edges, stretched point contacts, and intrinsically quasi-low dimensional compounds. 

1.1.3 Models in second quantization 

Let us now define the basic theoretical models for quantum wires. The Hamiltonian for an 
arbitrary band-structure of non-interacting fermions as shown in Fig. [7] can always be written 
in the formalism of second quantization 



system and the annihilation operator Cfc jCr annihilates a particle. These operators are defined 
for each spin index a=] ,j and each band separately and are summed over independently in 
the Hamiltonian. Spin and band indices will be suppressed in what follows until section 2.6. 
For simplicity we can assume periodic Born-von Karman boundary conditions over a finite 
length I = Na, which leads to discrete k values of k = 27rn/N. The usual anti-commutation 
relations for fermion operators must be obeyed 




(1) 



Here the operator c k creates a single electron \k,a) 



= cl 1 0) in a Bloch eigenstate of the 




(2) 
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Figure 7: Non-interacting electron dispersion. 



This ensures that the wave-function that is created by these operators is automatically 
anti-symmetric under particle exchange. Also the Pauli Exclusion Principle of no double 
occupancy c\,c\, = is obeyed. The number operator rik = c\ck in the Hamiltonian counts the 
number of fermions in the state \k) with eigenvalues of and 1, i.e. 



4 c fc 



0) = (3) 



and 



<kl*> = ck4io> = (4 {c fe ,4} - 44c fe )|o> = 4io> = \k) (4) 

The Fourier transforms of the Bloch operators correspond to field operators 

= 4E e ^ c *> ( 5 ) 



which also obey canonical anti-commutation relations. Here Xj=ja with j = 1,..,N labels 
the locations in the one-dimensional lattice with lattice spacing a. Therefore, ip^(xj) creates 
an electron in a localized Wannier state at lattice site Xj. The field operator is periodic 
4>{xj) = ip(xj + N) since e tkN = 1 for k = 2-irn/N . The inverse Fourier transform is given by 



Cfc 



1 N 

-J2e- ikj i>(xj). (6) 

3=1 



Example: Tight binding model 

In order to illustrate the use of second quantization in a simple model, we consider hopping 
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Figure 8: The interaction corresponds to a scattering process. 

between neighboring orbitals along a chain. The field operator ip{ x j) now creates an electron 
in the orbital wave- function (j){r — fj) , where fj is the three dimensional coordinate of the 
lattice site Xj. The so-called tight binding Hamiltonian is represented by overlap integrals 
t= — J d 3 rcj)(fj)AH(p*(fj + i) between neighboring orbitals, which allows transitions of the 
electrons ("hopping"). In second quantization this Hamiltonian is given by 

H = -tY J {^{x j )i>{x j+1 ) + tffa+xMxjj) (7) 

3 

Using the Fourier transform (|5|), we get 

H = "4 E E E (e- ikj 4e ik ' U+1) Ck> + e- M ^c{,e^c k ) 

j k k' 

= _ t Y j (e ik clc k + e- ik c\c k ) (8) 

k 

= —It n k cos k 
k 

where we have used the identity e%k ^ = which follows from the geometrical sum 

q 3 = — f^z — . This is an illustration of a band structure e k = —2tcosk with cosine depen- 

j=0 

dence where n k = c* k c k counts number of fermions in the state k. 
Interactions 

The understanding of electron-electron interaction effects is the main motivation for bosoniza- 
tion. The standard model of a generic density-density interaction is represented by a two 
particle operator. If the interaction potential between electrons in Wannier states at distance 
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Figure 9: Two tunneling examples: scanning tunneling microscopy and a conducting wire 
between two electrodes 



m is given by U(m), the corresponding interaction Hamiltonian can be expressed in second 
quantization using the fermion density operator p(xj) = ip^ (xj)ip(xj) as 

N N 

j = l 771=1 

N N 

= ^EEE4' e "" j/3C/ H E c{„c k „,e- ik "^e* k '"^ (9) 

j=lm=lk,k' k",k"' 

The sum over j results in a delta function with the condition k — k' = k'" — k" = Ak, so that, 

1 N 

ff « = jyE E clc k _ Ak U(m)c{,,c kl , +Ak e^ km 

m=l k,k",Ak 

= M E C k c k-Ak U ( Ak ) c k" c k"+Ak ( 10 ) 
k,k",Ak 

where we have used the Fourier transform and defined U(Ak) = ^ e l ^ km U (m) . The two 

777 

particle interaction therefore corresponds to a scattering process as represented in the diagram 
in Fig.[8j Particles Ck'-Ak an d c^-Ak ar e annihilated and after exchanging Ak two new particles 
ct, and c k are created. 

1.2 What is measured? 
1.2.1 Tunneling 

Tunneling is a common experimental setup for quantum wires. One example is scanning 
tunneling microscopy (STM), where tunneling occurs across the gap between the wire and the 
STM tip. A second example would be a ID wire, that is weakly connected to two leads at 
the ends, where tunneling may occur. In those experiments the current is determined by the 
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Figure 10: Schematic diagram of tunneling from a lead/tip to a wire 



tunneling rate, which is given by Fermi's Golden Rule 



r + {u,x) 



E' 



-pE„ 



n 



m 



5(lo - E n + E m ). 



(11) 



Here the probability of the transition from state \m) to \n) by adding one particle is con- 
sidered. The tunneling Hamiltonian between the lead/tip and the wire is assumed to be 
H = —t (ip* (x)i[>i ea ds + h.c.) and Z is the partition function. The tunneling current can be 
calculated from the tunneling rates by considering the empty and occupied states relative to 
the Fermi level given by the Fermi-Dirac distribution, f(uS) = 1/(1 + e^ w ) 



l(y,x,0) = el dujp leads (uj-eV)[f(u J -eV)T + (u J ) - (1 - f (u - eV)) T~ (u)] 

(Tunneling in) — (Tunneling out) (12) 



as illustrated in diagram [TUJ It is useful to express the tunneling rate in terms of the so-called 
local spectral weight, defined as 



l + 



-PEr, 



n 



2 1 + e _/3w 

5{u-E n + E m )= ?— 2 r+, (13) 



assuming finite temperatures (5 = 1/ksT. Inserting (|13f) into the expression for the current 
(|12p and performing straightforward calculations using the Fermi-Dirac distribution, we arrive 
at a simple equation for I(V,x) in terms of the local spectral weight 



I(y,x,p) =27rt 2 e / dwp leads {u-eV){f{uj-eV)- f{u))A{u,x). 



(14) 
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Finally, in order to calculate the local spectral weight A(u),x), we can rewrite the delta 
function as 5(u) = ± £™ e lut dt = Re \ J °° e lwt dt, so that from (fT3|) 



A(u,x) = ±Y,( e ~ f}En + 



-PEr, 



n 



ijr (x) 



m 



2 1 r°° 
R e - e i(^-E n +E m )t 

7T J0 



dt (15) 



Then, using e tHt \m) = e tEmt \m), we arrive at the Green's function representation 



iE m t 



A(u,x) =Re ( -^Yl [°° dteiuJt (e~ pEm + e^ En ^ (n $ (m \e iHt ^e- iUI 1 

\ n,m 
1 f°° I 

= Re- / dte wi (V(x, tW(x,Q) +ipUx, 0U(x,t) 

7T Jo X 



n) 



1 



Im / G H (t,x)e luJt dt 
n Jo 



(16) 



where G R (t,x) is the retarded Green's function 

G R (t,x) = -i({^(x,t),^(x,0)})e(t). (17) 
Here the thermal expectation value (^4) = ^ e -PE n ^ is used. At zero temperature 



the local spectral weight reduces to the local density of states (LDOS) 

p(u,x) = V Urn ipUx) 0) S(u - e m ) = -Im / G R (t, x)e iu}t dt 
™ 1 -/o 



and the current is 



eV 



I(V, x, (3) =2irt e du; pi ea ds{u ~ eV)p(u, x) 
Jo 



(18) 



(19) 



Since Green's functions such as (|17p play a central role in condensed matter physics, it will 
be one of the main goals in the many-body theory to calculate time correlation functions. 



1.2.2 Photoemission spectroscopy 

Photoemission spectroscopy (PES) is a standard technique to determine the electronic proper- 
ties of condensed matter systems. Photons with definite energy and direction strike an object 
and electrons are emitted via the photoelectric effect as shown in Fig. QTJ In return, the 
intensity of the electrons as a function of absorbed energy and momentum gives information 
about the object. In inverse photoemission electrons are used to probe the sample and the 
emitted photons are analyzed correspondingly. Interestingly the mathematical description is 
very similar to the tunneling processes considered above. Just as before, we can argue that 
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wire 

Figure 11: Schematic diagram of photoemission 

the emission rate of photons T + in inverse photoemission at a given energy and wave-vector 
is related to the probability for adding a corresponding electron into the system. According 
to Fermi's Golden Rule 

r+(u,q,/3) oc i^e- /3 ^|(n| C t|m)| 2 5(u;- J E; n +S m ), (20) 

n,m 

where now the probability of the transition from state | m) to | n) by adding one particle with 
wave-vector q is the relevant quantity. In this case it is useful to define the angle resolved 
spectral density 

2 5(uj-E n +E m ) <xT + (uJ,q). (21) 

n,m 

The angle integrated spectral density and the space integrated local spectral weight in 
equation (|13p are always the same and represent the total spectral weight as a function of 
lo. Following the analogous calculations as in the previous section we can calculate the angle 
resolved spectral density in terms of a retarded Green's function in momentum space 

1 f'°° 

A(u,q,p) = Im / dte lult G R (t,q,f3), (22) 

where 

G R (t,q,P) = -i({c q (t),cl(0)}). (23) 

1.2.3 Conductivity 

When speaking about the physics of "wires", it is natural to think about the corresponding 
conductivity. However, in order to apply a voltage and measure the current, perfect contacts 
to external leads are required, which cannot be realized in most experiments with truly one- 
dimensional wires. Nonetheless, let us imagine a setup in which electrons are adiabatically 
guided into a single-channel wire of length £ without any scattering under an applied voltage 



A(u,q,/3) 



i + e 



-PEr, 



n 



m 
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Figure 12: Schematic diagram of metallic leads and a ID wire 



V as shown in Fig. [T2l 

The standard theory of linear response gives the general conductivity in terms of the Kubo 
formula 

a(q,u) = - / e^([J(q,t),J(-q,t)])- (24) 

Naively the DC conductivity is then simply given by cr(0, 0), but since we are dealing with 
a finite wire it turns out that it makes a difference if we first take the limit to — > (current 
oscillating between the leads) or the limit q — > (current oscillating within the wire). In the 
latter case we obtain 



-i pt poo 

G(w) = Mma(q,u) = — dx dte^ 1 ([J (x,t), J '(0,0)1) . (25) 

For interacting systems the outcome actually depends on the order of limits. However, we 
will not address this question here. Instead we consider a perfect ballistic non-interacting wire 
in order to show that even in this simplest case the conductivity is finite and quantized. In fact, 
we can derive the conductivity without the calculation of correlation functions by considering 
the semi-classical acceleration of electrons in the field of an applied voltage E = V/l. The 
crystal momentum changes with the applied force hk = —F = eE = eV/l. The current is 
given by the number of electrons in current carrying states per unit time. The number of 
electrons that contribute to the current is simply given by counting the states in an interval 
Ak as # of e~ = Ak£/2ir (see Fig. [T3"|h Since electrons are constantly ejected into the 
right lead, the number of electrons in the current carrying states remains finite and can be 
determined by the number of electrons that are accelerated into those states per unit time via 
the following simple calculation 

#o/ e~ Ak i ■ I e 2 V 

I = e 7 = e-r— — = ek— = ——r = GV 

At At 2ir 2vr 2irh 

where G = e 2 /h is the quantized conductance per ballistic channel (2e 2 /h for electrons with 
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Figure 13: Electrons are accelerated by the voltage into current carrying states. 

spin) . The corresponding energy is dissipated in the right lead where the accelerated electrons 
quickly reach thermal equilibrium, while the wire remains cold. In general it is also possible 
to take a finite transmission and several channels into account, which gives 



where < t n < 1 are the eigenvalues of the transmission matrix T^T between the n different 
channels. This is the famous Landauer-Biittiker formula for quantized conductance of wires 
and point contacts. 

2 Bosonic description 

2.1 Linearization of the fermion dispersion 

A bosonic description of the one-dimensional wire will be useful in the treatment of interac- 
tions. Let us first show in detail how a non-interacting fermionic system given in equation ([1]) 
can be expressed in terms of a bosonic Hamiltonian in second quantization. For this equiv- 
alence to work we require the dispersion relation of the fermions to be linear in the range of 
interest. This is always approximately true for excitations with small energy around the Fermi 
points. Therefore, we will restrict ourselves to consider low enough energies so that we can 
describe the band structure in equation (pQ) by an energy that depends linearly on the wave 
vector k, which is the only approximation we will make. The starting point is an arbitrary 
free electron dispersion in equation ^Q) 




2e 2 



channels n 




(26) 
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Figure 14: Allowed k-values around ±kp 



At zero temperature all states between the Fermi wave-vectors -kp to k F are occupied. The 
number of particles in the Fermi sea is therefore given by 

n = k F N/ir - 1 (27) 

if no is even, or no=k F N/tt if no is odd (times two for electrons with spin). 

Since we are only interested in low energy excitations, we now restrict the allowed wave- 
numbers to a range A around the Fermi points — k F — A < k < —kp+A and kp— A < k < kp+A 
as shown in Fig. Q3] corresponding to right- (kp) and left-moving (-kp) fermions, respectively. 
The range A is determined by the region over which the band structure is approximately linear 
and depends on the model (typically corresponding to an energy range of about 1/10 of the 
band width). Keeping only states in this range the linearized effective Hamiltonian is given 
by 

kp+A -k F +A 

H PS ^ e kc[ck + ^2 e k c \°k ( 28 ) 
k=kp—A k=—kp—A 

In this range the energy dispersion can be expanded around ±kp 

for k ps k F (29) 

- " Kp 

for the "right movers" and analogously 

for k ps -kp (30) 

" " — Kp 

for the "left movers". Here the Fermi energy £-k F = £k F = must be zero since all occupied 
states have negative energy and all empty states have positive energy at zero temperature. 



+ o 



(k-k F y 



e k ps e-kp + (k + k F ) — ^ 



+ o 



{k + k F f 
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Defining a Fermi velocity 



vf 



dk 



kp 



' dk 



(31) 



-kp 



we can therefore write e k ~ vp (k — kp) for right movers and e k ~ — vf (k + kp) for left 
movers. It is useful to define new quantum numbers \k\ < A shifted relative to the Fermi 
points ±kp and corresponding left and right moving fermion operators 
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Ckp+k 
c -k F +k 



(32) 



In this way the effective Hamiltonian has the simple form 

A 



£ vpk (cf 4 - c ^cl) 



(33) 



k=-A 



Accordingly, the original fermion field operator in equation also splits into a left and a 
right-moving part 



-kp-A 



N 

K 

1 / kp+A 
= ^ \kS-A 

= ~^Y1 (e lkFj e ikj c§ + e- ikFj e ikj 4 



e ikj c k 



+ £ < "'7. 



k =A—kf 



(34) 



) 



where we have defined left and right moving fermion fields on the length of the wire £ = Na 



To S C k 



fl ikxj/a 



V£ 



k=—oo 
oo 



L ikxj/a 



c k e 



(35) 



k=— oo 



In those definitions it is useful to extend the range of summation to infinity, because in this 
way an inverse Fourier transform can be defined as a continuous integral 



i' R = ^f Q e- ikx '^ L/R {x)d. 



x 



(36) 
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One may object that we have now changed the range of summation to infinity in equations 
(I34p and (|35p . since the cutoff A is supposed to be much smaller than the bandwidth. However, 
remember that A also represents an upper limit for the validity of the effective low energy 
description. Indeed, extending the states of the linear dispersion relation to infinity or leaving 
them out makes no difference if we always restrict ourselves to low energies, since those states 
will then never take part in any physical excitation. Therefore, all states with fc>A and -k<A 
are unphysical anyway. By taking the summation range to infinity we have effectively taken a 
continuous limit for the field operators, which is mathematically more convenient. However, 
this does not correspond to an additional approximation. Note that the anti-commutator 

{^(x)^M2/)} = |E e " jWvfcVa { c f'^} = 7 £ e-^-y^ = 5(x-y) (37) 

kk' n=—oo 

is now normalized as a delta function. In general, the normalization would actually depend 
on the choice of the cut-off A, but by including the non-physical states a more convenient field 
operator with canonical anti-commutation relations has been defined. 

2.2 Excitation spectrum 

In order to understand how the fermionic spectrum can be represented by bosons we will 
consider individual excited states as e.g. depicted in Fig. [15] on the right moving branch. In 
terms of fermionic creation and annihilation operators the state is expressed as (omitting the 
index R for right movers temporarily) 

cA4cU4\GS) , (38) 

where we use integer labeling relative to the highest occupied state of the Fermi sea for 
simplicitj0 

n = kN/2ir + 1/2. (39) 

Note that the order of the operators matters, since they anti-commute. Alternatively, we can 
construct this state in two parts: 

1. First we create four additional fermions in the lowest unoccupied states c\c\c\c\ \GS) 
as shown in Fig. [TBJ The energy cost for this is given by 1 

2ttv f (1 3 5 7\ _ 2ttv f 16 

1 Note that the Fermi level can be chosen to be in the middle between two states if no is odd. The highest 
occupied state n = is at fci? — n/N and k is measured relative to kF- 
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Figure 15: Left: The ground state of the filled Fermi sea \GS). Right: Typical excited state. 
The dots represent occupied states, crosses unoccupied states. 

or in general by 

En R = ^fn%, (40) 

where ur is the total number of additional right movers. The energy cost is the same for 
removing —tir particles. 

2. Secondly we create particle-hole excitations by shifting up the individual fermions 
(preserving their relative order). In our case we have to 

• shift the top fermion up by 4 steps with energy cost 1 x 4 x 2ttvf/N , 

• shift zero fermions by 3 steps, 

• shift the two next fermions by 2 steps each with energy 2 x 2 x 2ttvf/N 

• finally shift the two next fermions up by 1 steps with energy 2 x 1 x 2itvf/N 

By this shifting procedure any particle-hole excited state can be created. If we now interpret 
that each shifting by n steps corresponds to a boson of level n, the particle-hole excitations in 
the part 2 above would correspond to 

• one excited boson on level 4, with energy 1 x 4 x 2ttvf/N 

• no bosons on level 3 

• two bosons on level 2 with energy 2 x 2 x 2itvf/N 
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Figure 16: Excited state after adding four fermions in the lowest unoccupied levels 
c\c\c\c\\GS) after step 1. 

• two bosons on level 1 with energy 2 x 1 x 2-kvf/N 

In this way any particle-hole excitation can in principle be represented by bosonic occupation 
numbers. This hand-waving picture is almost the entire secret to bosonization. The full story 
involves only a slightly more complicated linear combination of shifted states as we will see 
now. 

2.3 Bosonic operators 

From the previous section it is intuitively clear that "fermion shifting operators" may be 
represented by bosons. In this section we will make this mathematically precise. A good 
definition for a fermion shifting operator is given by 



This operator changes the wave- vector of right moving fermions by an amount k, summed over 
all possible states k\ resulting in a superposition of shifted states. This is a slight difference 
from the intuitive operations discussed in previous section where we have only shifted the 
uppermost fermions at a time (part 2), but clearly equation (|4ip amounts to the same basic 
idea. 




(41) 



k> 
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The operator is not hermitian. In fact, the hermitian conjugate changes the /c-vector by 
the opposite amount 



+00 

W?) f = £ c v4+k = Y, c k>-A = p-h (42) 

£:'=— 00 k' 

If p R with k > creates an excitation as suggested in the previous section, the hermitian 
conjugate annihilates it, which is indeed what we would expect of bosonic operators. In the 
last equation we have shifted the summation variables, which is unproblematic because we 
have used the convenient trick of extending the summation over an infinite range as discussed 
above. 

The final and most important step in order to relate the shifting operators in (I4ip to bosons 
are the commutation relations. Using the definition (|4ip we get 



(43) 



\ n R n R] _ V (r R ^ r R r R ^ r R r R] r R r R] r R \ 
[P-kiPk'l - \ C k"-k C k" c k"'+k' C k"' ~ c k"'+k' C k"' c k"-k C k" J 

k",k"> 

- V (r R] ir R r flt \ r R r R ^ { r R r R] \ r R 
~ 2^ \ C k"-k \ C k"i c k"'+k> j C k"' ~ c k"'+k' \ C fc'"> c k"-k / C k" 

k",k"' 

_ «t R] r R r R _j_ R] R] r R r R \ 

c k"-k c k"'+k' c k" L k"' ' L k'"+k' c k"-k c k'" c k" J 

- V (r RJf r R r R] r R \ 
~ 2^1 \ C k"-k C k"+k' ~ c k"-k+k' C k") 

k" 

= 0, 

where the two terms on the third line cancel and we have again shifted the summation variable 
in the last step (k" — ► k" — k). Therefore we find that the shifting operators generally commute. 
However, if k = k' , the summation variable cannot simply be shifted because in that special 
case the result corresponds to the subtraction of two infinities 

[p R k, Pk } = £ ( nfc "~ fc ~ nk ") = 00 ~ 00 ( 44 ) 
k" 

Indeed any arbitrary result can be produced by subtracting two infinities! However, we re- 
member that these are not real infinities (coming from non-physical states below the lower 
cutoff k<-A ). Therefore, using the fact that all states below the lower cutoff always have to 
be occupied, we can cancel the unphysical infinities by using the cutoff at —A as a reference 
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point, since we require nk = 1 for k < — A 

[p-k,Pk] = ^2 ( n k"-k-n k ")+ ^2 (, n k"-k-n k >>) 

k"<-A k">-A 

= ^2 nk " ~ ^2 Uk " 

k.">-A-k fc">-A ^g^ 

= ^ n k " 

-A>k">-A-k 

_ kN 
~ ~2T 

where all terms in the first sum on the first line are zero and we have used the finite level 
spacing k=2'Kn/N. In summary, we have 

[p R k,P§]=hk>^ (46) 

which indeed corresponds to bosonic commutation relations. 

For the bosonic operators to be useful we still must relate them to the Hamiltonian. Using 

MOOjl ii — upf 

k' 

given by 



equation ([33]) H = Y2 v F k ' [c k ,c R — c k Jc k ,j, the commutator with the shifting operator is 

k' ' 



R\ R R] R 
C k> C k' ' C k"+k C k" 



[H,p§] = 

k',k" 

V wk' (r R Ur R r R] \ r R r R] { r R r R] \ r R 
- V F K \ c k> \ c k' ' c k"+kf c k" ~ c k"+k \ c k"i c k' ] c k' 

k',k" 

= ^2 v Fk' (c k )c R _ k - c k } +k c k , ) 

k' 

= ^2 y p U k ' + k ) c kKk4 - k ' c k>+k c k' 



k> 

= v F kp k (47) 

It turns out that the two equations (|46p and (|47p are in fact sufficient to determine the 
boson algebra and the bosonic representation of the Hamiltonian. In order to see this better 
we define conventional creation boson operators for positive fc-values (i.e. shifting particles up 
k > 0) 
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while the hermitian conjugate annihilates an excitation by shifting particles down 



For left movers negative fc-values correspond to an excitation, so that for k > 

p L - k = E c tLA> = *Mb? (so) 

k' 

The phase in the definition is in fact arbitrary, but we choose ±i for later convenience. Inserting 
this definition in (|46p and (|47p we recognize the canonical commutation relations 

= $kk', (51) 
= v F kbf (52) 

and likewise for left-movers. As shown in the treatment of the quantum harmonic oscillator in 
any quantum mechanics book, operators and b with such commutation relations are called 
creation and annihilation operators and can be used to build up the entire spectrum. In our 
case, we have such an oscillator spectrum for each k > separately, describing the particle-hole 
excitations. Now, also taking the energy for adding/removing extra particles from equation 
(|40p into account (see Fig. I16p . we arrive at the complete Hamiltonian in bosonic form 

H = v F £ k (bP% + tfttf) + ^ (n| + n\) (53) 

k>0 

where (n 2 R + n^) nvp/N is the energy in the first step of adding particles expressed in count- 
ing operators tir and m (so-called zero modes) and vpk (bj^bj* + b^bj? J corresponds to the 
second step of shifting of particles (so-called oscillator modes). 

2.4 Fermion operators in terms of bosons 

We have achieved the most important step of expressing the Hamilton operator entirely in 
terms of bosons. However, if we wish to calculate any physical expectation value we must 
also be able to express more general fermion operators in terms of the bosons. The bosons 
are well defined in terms of the fermions by equations (l48H50p . but we now seek the inverse 
transformation. 



°k > °k' 



H,bf 
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2.4.1 Left- and right-moving fermion densities 

Using the concept of left-and right moving fermion field operators in equation (|35p we can 
immediately define a corresponding left and right moving fermion density 



p R ( x ) = ^ R ( x )^ R (x). 



(54) 



Technically this operator is divergent because the ground state expectation value is 



x/a 



(<«) = \ ± (cf 4 




= oo 



(55) 



k,k' 



fc=— oo 



where we have used (|35p . However, we understand again that this infinity comes from intro- 
ducing unphysical states below the lower cutoff k < —A as discussed above. Since those states 
are always occupied, the same divergence will appear whenever this operator is applied on 
any physical state. However, we can simply remove the divergence by subtracting the ground 
state expectation value 



This procedure of subtracting the ground state expectation value is called normal ordering 
and is indicated by the dots at the sides. This definition is indeed useful because we are 
actually only interested in excitations on top of the well-defined ground state, namely the filled 
Fermi-sea. The normal ordered left- and right-moving densities therefore measure density- 
fluctuations relative to the ground state. Instead of subtracting the ground state expectation 
value, we can equivalently always re-order the plane- wave fermionic operators and c k , so 
that the annihilation operators are to the right for k > and to the left for k < while keeping 
possible minus signs from the anti-commutator, i.e. 



Note, that the normal ordering can be omitted if k ^ k', since : c k cjp : = c k c k , = —c k ,c k 
already follows from the regular anti-commutation relations in that case. For normal ordered 
left movers the annihilation operators are to the left for k > and to the right otherwise. 
We now consider the Fourier transformation of the left- and right-moving densities, first 




(56) 




for k > 



for k < 



(57) 
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assuming k ^ 0. By using (|35|) and J Q dx e %kx l a = £5k t o we find 

e ikx / a dx : p R (x) : = - f dxJ2 
J Q k"k' 



ikx/a —ik" x / a ik'x/a . R\ R . 



/ > • C k+k' C k' ■- C k+k' C k' 



(58) 



On the second line we used the fact that the normal ordering in (j57[) does not affect the 
expression for k ^ 0. We therefore find that the Fourier components of the right moving 
densities correspond exactly to the fermion shifting operators defined in (|4ip ! 
For k = we have 



Jo 



dx 



: p R {x) : = \ fdx ^ e-WW* : c$cg : 
- 2_^- C k' C k' ■ 



(59) 



/ , C k> C k' ~ 2^ C k' C k' 



k'>0 
n,R 



k'<0 



We therefore find that the zero mode : p R : corresponds to the total number of right moving 
fermions relative to the ground state : p R := n R (the third line corresponds to the number of 
excited particles k'> minus the number of holes k'< 0). In terms of the bosons in (|48p the 
equations ([351) and ([39]) can be summarized as 



[ i dxe ik ^:p R (x): = Y,--c B J +k c§:=< 
J v 



kN . 
2tT 



nR 



b R] for k > 



for k = 



(60) 



kN, 
2~7 



n/^-bf forfc<0 



We now make use of the inverse Fourier transform of equations (|58[) and (|59p . which is given 
by (fc = 2vrn/iV) 



: PM^J : 



= 7 E : Pk ■ = ]Y,\M ( ib * eikX/a ~ ^ e ~ lkX/a ) + T 

fe fc>0 

1 oo 



(61) 



n { ib R e iJ r nx - ib^e 



71=1 
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For the left moving density we can make the analogous calculations to find 




In summary, we have derived central bosonization formulas which can be used in order to 
express fermionic densities in terms of bosons. 



2.4.2 Boson field operator 

It turns out to be useful to introduce a boson field operator <j>R in order to express the 
bosonization formulas (|6ip and (|62p more compactly and also in order to make contact with 
conventional bosonic field theories. Let us the define the fields 4>r and cj)L 



OO -. 



n=l 



in terms of the boson operators bn and b^ L as defined in equations (|4"8|) and (f50l) . The 
operators = ^pnn R and Ql = y/^fiL measure the right- and left-moving particle number 
tir and til- The operators 4>q and 0q are defined as the canonical conjugate to the number 
operators 

[&Qr] = ~\ 

[<^,Ql] = \ (64) 

The definition in equation (|63p is useful, because we can immediately express the fermion 
density in equations (fBTj) and (|6"2j) as the derivative of the boson field 4>r and c/>l 

■ Pr{x) : = -^d x (p R (x) 

f (65) 

: Pl{ x ) ■ = —^d x (j)L(x) 
\/7T 
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Moreover, the Hamiltonian can also be expressed in terms of the boson field as 
Inserting the definition and using j L dxe l ~ nx = S n< o we find that 



H = av F dx {{d x 4>R? + (dM 2 ) (66) 
Jo 



H = 2 ^fY. n ( W + W) + v -§(Ql + Ql) (67) 

n=l 

which indeed is identical to the boson expression we found in equation ([53]) . 
2.4.3 Fermion field 

In order to express arbitrary fermion operators and calculate Green's functions, it is necessary 
to describe the fermion fields ip R (x) an d V'z,( a; ) m terms of bosons. We know that a right- 
moving fermion field ip R (x) annihilates a fermion locally at position x and therefore changes 
the density Pr{x) = ip R (x)ijj R (x) at that point. This action is reflected by the following 
commutation relation 

[ip R (x),p R (x')] = ip R (x)ip R (x')ip R (x') - ip R (x')tl) R (x f )i/j R (x) 
= [^ R {x)^ R {x')]^ R {x') 

= 5{x-x')^ R {x) (68) 

Our goal is therefore to find a bosonic expression for ip R {x) which reproduces the commutation 
relation f)68f) . The bosonic expression for the density : p R (x) : is already known from equation 

AMD- 

Before we go further, let us review a general relation for commutators of two operators A 
and B. Assuming that [[A, B] , A] = we can show that 

OO -. OO -. 

B ] = E - B ^ = E ~ nAn ~ X l A > B ^ = eA l A > B ^ ( 69 ) 

n=l n=l 

Applying this equation to typical boson operators, i.e. A=b, B=b + with [A, B] = 1, we get 



ae Qb 



(70) 



The commutator with an exponential of boson operators is again proportional to the exponen- 
tial which is approximately what we want in equation (|68p . We are therefore motivated to try 
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the following ansatz to describe the fermion field with an exponential of a linear combination 
of boson operators with arbitrary coefficients a n and f3 n 



Using 



b R b m 



ip R (x) oc exp I [ a n(x)b R + Pn{x)b R 

\n>0 

5 n y and for each index n,n\ we get 

ipR(x),b^ = a n ip R (x) 
[i/) R (x),b%] = -Mr(x) 

Together with equation (|6ip we can therefore calculate the commutator 



(71) 



(72) 



ty R {x),p R {x')] =iY J ^{e^ x ' [^ R (x)X]-e- l ^ x ' 

n=l 
1 X 



(73) 



n=l 



So [V>i?(x), Ar(x)] is indeed proportional to Tp R (x) as desired. According to ([68]) we must choose 
the coefficients so that we get a delta function 



(/3ne i2 <" x ' + a n e 'T* 1 ) = £(x - x) 



(74) 



n=l 

+oo 27r 

Comparing with the equation | ^ e -«— «(^~^ ) _ j^, _ ^/^ wg conc i U( ^ e that 

n=— oo 



(75) 



This reproduces the delta function almost perfectly, up to the n = term. We therefore guess 
that we need the zero mode in the exponential ansatz (|7ip as well. Prom equations (|64p and 
(1691) we know 



Therefore, 



oc exp ( iV^<p R + ^ A={e.- i2 f nx bf + e i2 r nx b R ) 

V n>0 ^ n / 



(76) 



(77) 
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fulfills the relation 



[ip R (x),p R (x)] 



Qr 



+ \ 0^*' M*), 6ffl - ) 

n=l 



71=1 / 

= 5(x - x') iI>r{x) 

(78) 

as required in equation (|68|) . Including the operator <5i? i n the exponential (|77"|) does not 
change this relation, so that we can arrive at a compact form by comparison with (|63p 



ipn(x) oc exp ^i\/47r0i?(x)^ (79) 

which is the famous bosonization formula for fermionic fields. Following the analog calculations 
for left-movers we find that 

V'zOc) oc exp ^— za/Itt^l^)) . (80) 

These two formulas will be useful when calculating Green's functions and general correla- 
tion functions as we will see later. 

2.4.4 Field commutators 

So far we have not specified the overall proportionality constant of the bosonization formulas 
(|79p and (|80p . This can be fixed by using the proper normalization introduced in equation 
(|37p . but it must be noted here that even that normalization can in principle be chosen to be 
cutoff dependent, so that we could equally well leave the overall constant arbitrary, which is 
often done in the literature. Nonetheless, calculating the anti-commutator in bosonized form 
and calling the proportionality constant C, we get 

\ipUx),ip R (y)) = C 2 (e-^^le 1 ^^} 

V- / x ( 81 ) 

= C 2 e -iV4^{<j} R {x)-4> R {y)) ( e 2-K[<t> R (x),<t> R (y)} + e -27T[^ R (x)^ R (y)]^ 



where we have used the so-called Baker Hausdorf formula 

e A e B = e A+B e \[A,B] = e B e A e [A,B] (82) 
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(assuming [[A, B] , B] = [[A, B] , A] = ). Hence we need to calculate the bosonic commutator 
(using flggp and ([MD) 



m ( x ~y 

i 



+ E 



2^ 



7^_ 1 47rv nn' 

oo 

+ E 



e i^(ni-n't/) 



b R b Rt 



_|_ e -i^f{nx-n'y) 



n=l 



47rn 



= ^+-Lln(l-e-^^))-^ln(l 
2i 4vr V 7 4vr V 

where we have used the expansion of the logarithm 



Since 



hi 



E 

n=l 



l-e~ 
1 - e l 



n 



we have 



= -ln(l-e- a ). 

i (it — a) for a > 
for a = 
_ i (— 7T — a) for a < 



- for x — v > 
4 y 



*>R (x) , 4>r (y)\ = - sign(x - y) = < 



for x — y = 
— — for x — y < 



Inserting this in equation (|8ip we get 

{^(x),^(y)} 



for x / y 
2C 2 for x = y 



(83) 



(84) 



(85) 



(86) 



This expression agrees with (j37|) if the proportionality constant C 2 = S(0)/2 is chosen to be 
a delta function infinity (or alternatively is cutoff dependent, depending on the summation 
range in (l37l) ). 

Following the same steps for left movers, we find analogously 



[<f>L (x) , <t>L (y)} = sign(x - y) 



(87) 



Similarly, the anti-commutation relations between right- and left movers can be guaranteed 
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by demanding suitable commutation relations between the zero mode creation operators 







which yields [<fi R (x),<fi L (y)] = | and therefore {i()r(x), 4>l(i/)} = according to equations (f63l) 
and (|79p . Alternatively, the commutation relations between left- and right-movers can also be 
fixed by using so-called Klein factors, which can also handle several fermion channels. 

2.4.5 Bosonic excitation and zero modes 

At this point it is useful to pause and look at the excitations of the bosonic Hilbert space 
in more detail. For the bosonic "oscillator" modes it is clear that the ground states of the 
Hamiltonian in (|67p simply corresponds to the vacuum |0) which can be excited with an 
arbitrary number of bosons with the creation operators b n and b n for each n separately. The 
vacuum is defined as 

b n \0) = for all n (88) 

The total fermion number relative to the ground state can be added and removed with the 
zero mode operators exp (— iy/^jrcj)^ and exp u'V47r<^) , respectively. The corresponding 
relation n^e -1 ^^^ |A) = e ~ ly ^<t>o (n R + 1) |A) for any state |A) can be shown with the help 
of equation (|76p and is left as an exercise here (remember that ur = Qr/\/tt is the counting 
operator). The fermion number creation always acts from below, i.e. the entire corresponding 
state is "pushed up" one step by acting with exp (— i\/47r0^) including any possible particle- 
hole excitations. An arbitrary number of fermions can be removed or added, but only in 
integer numbers (see also Fig. [To 7 ]) . The mathematical reason for this is actually the periodic 
boundary condition tp(xj) = ip(xj + ]y) . According to equations (|34"]) and ([35]) the boundary 
condition implies for the left and right moving fields 

1>r{x) = i> R {x + e)e ikpN , i> L (x) = ip L (x + i)e~ ikFN (89) 



This translates into conditions on the boson field (|63|) by the use of the bosonization formula 
(|79p . Since the sum of oscillator modes in the boson field operator ^ 1 (e l ~ x b^ + 



n=l 

_„-2irn j Rf 



e * e x b n ) is already periodic in £, we are left with the following condition on the zero modes 
(using (|82p since 4>q and Qr do not commute): 

> I — e iV^QR e -i^ e ik F N 

— > V^ttQr = 2vrn - k F N + vr (90) 
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We see that the number operator Qr/^ = ur must indeed be an integer, since kpN = 
(no + l)vr is an odd multiple of ir according to equation (|27p. Therefore, zero mode parti- 
cle excitations exp (— i\/47r</>Q^) can only be added/removed in integral numbers. Arbitrary 
real numbers in the exponent would lead to unphysical states outside the Hilbert space. In- 
terestingly, the spectrum and Hilbert-space of the zero modes therefore have a one-to-one 
correspondence to a single particle on a ring, where Qr plays the role of the momentum oper- 
ator and 0o i s the corresponding position operator with canonical commutation relations as 
defined in (pi) . 

2.4.6 Summary of the bosonization formulas 

Before we continue, let us summarize the most important bosonization formulas. 



Linearization 



c k F +k 



oo 



fc=— oo 



fc=— oo 



Bosons in terms of fermions 

~iM h k =Pk= Jo : ■ etkx/adx = E 4lk4 ^ * > 



= & = Jo : pfc) ■ e ~ lkx/adx = E HA 

v k' 
Qr/l/V^ = riR/L = Jo ■ Pr/l{x) : dx 
H = v F Y,k (fif 6« + 6jft 6 £ ) + ^ (n 2 + n | } 

fc>0 v 7 



Fermions in terms of bosons fields 

: Pr(x) : = -^d x 4> R (x) : p L {x) : = -^d x (f> L (x) 

4>r{x) oc exp {i\f^n(j)R{x)) iPl{x) oc exp {—i\f^K(j)i J {x)) 

H = av F J*dx {(dM 2 + (dM 2 ) 

where 

oo 1 



for k > 



n=l 
oo 



n=l 



47rn 



(USD 
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(® 
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2.5 Correlation functions 



As a typical application for bosonization, we would now like to calculate the correlation func- 
tion for the fermion fields. 

2.5.1 In space 

Correlations along the wire decay according to the following expression 



y (x)V(y)) ~ ae~ ik ^ x -yV a (f R {x)^ R [y)) + ae^-v)/* (ar)^(y)) (91) 
in the linearized approximation (|34p . The left and right movers are uncorrelated 



because of the zero modes: If a right-mover is added exp (— iyAftcf)^} it must also be removed 
again with exp (iV^r^^) in order to get the ground state back. Therefore, in particular 
(0|exp (-iv/47r^) |0) = . 

Using the bosonization formula in real space we can write for the right moving correlation 
function (also using ([52]) and ([85])) 

/ Vt? 0*0 1/>R (y)) =C 2 ( e -^(^(^)-^(s/))\ e 2n& R (xU R (y)] 

X ' ) r- [ (92) 

_ (j2 / e -iV4n((t> R (x)-<f> R (y)) \ e iirsign(x-y)/2 

In order to calculate the expectation value of the exponential, we have to consider the zero 
modes and the oscillator modes separately. The expectation values of the zero modes are 
simple since the addition/removal operators cancel in the sum 4>r(x) — (j>R.(y) and we are left 
with 

(0| exp (-iVI^Q R (x - y)/e) |0> = 1 (93) 



exp(a& + 0) o) = (0| e ab ^ e ph |0> = exp(a/3/2) (94) 



which follows from the fact that the ground state has no additional particles Qr |0) = 0. 

For the ground state expectation value of the oscillator modes we can use the Baker- 
Campbell-Hausdorff formula (|82p in order to arrive at the general expression 



where we have also used equation ([88]) . In fact the relation ([M]) is just a special case for T = 
of the more general cummulant theorem for bosons (exp(/)) = exp (/ 2 ) /2 (see appendix 
equation (j!33p ). which is valid for any temperature and any linear combination of bosons / as 
shown in the appendix. With the help of the cummulant formula (j!33p it is also possible to 



32 



calculate correlation functions at finite temperature and finite system sizes, which is explained 
in detail in Ref. [15] (also for "open" boundary conditions). For simplicity, we will restrict 
ourselves to the ground state correlation at T = 0. 

Equation ([92]) can be evaluated by inserting the mode expansion ([63]) into $4 



(^ R {x)ij R {y))=iC 2 s\gu(x-y)U A" J \ (95) 

Using again the expansion of the logarithm (|84|) . we find for the sum 



n>0 



= hi(l-e iS r^- y ^ +\n(l-e- i2 -r {x - y) ^ -2 lim (Ins) 
where the last term corresponds again to an infinity 



oo oo 

1 ,. \ - e 



E- = lim = — lim In (l — e £ ) = — lim In e. (96) 

n=l n=l 



However, this infinity combines with the prefactor C 2 , so that a finite result is obtained using 
the relation eC 2 =l/£. 



i>R ( x ) i>R (v)) = lim iC 2 sign(x - y) exp ( - In |- _ — 

/ £— >-o y i g J f i 3 ' y) \ 

i s\gT\(x — y) 
= 27|sinf(x-y)| 



£ 



(97) 



2^sinf(x-y) 



The relation eC 2 =l/£ can be shown, e.g. by evaluating the correlation function (|92p directly 
from the definition (|35|) by using the geometrical sum (and the quantization of the fc-values 
relative to kp as in equation (|39p ). The calculation for the left moving correlation function 
follows the same steps, except for a different sign in equation (IS] 



Therefore, in summary from equation ([91 



, + / s , / \ \ sin kp[x — y) . , 

V f x v (y ) = „ . ^7 y \ 99 

/ iVsm|(x-j/) 
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It is left as an exercise for the reader to show that this is indeed the exact correlation function 
which can also be obtained without any approximations directly from ([5]). 



2.5.2 In time 

The time correlation function plays an important role, e.g. in order to evaluate the Green's 
function (|17|) . Before we can evaluate the expectation values we have to determine how the 
operators evolve in time. Using the general formula 

A(t) = e iHt Ae- iHt (100) 

it is straightforward to show that 

Ml— ^-t 



2tyv p 



Qr/l^) = Qr/l 

4 /L {t) =^ L T Q R/L V -£t 

Inserting these results into equation (|63H we can immediately obtain a time-dependent mode 
expansion 



n=l » 
n=l 



(102) 



This is a remarkable result in many ways. First of all we see that the right- and left-movers 
are only functions of the right and left-moving light-cone coordinates x±avpt. Secondly, with 
the help of this result the calculation of any time-space correlation function is no more or 
less difficult than the calculation of a pure space correlation function by substituting the time 
dependent light cone coordinates for x — y in (|97p . That means that we have in principle the 
tools to solve any dynamic problem. In particular, we can immediately generalize the result 
(|97p to express the time correlation function 



where it must be remembered that t = lim (t — ie) in order to ensure that the summation 
converges. 
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2.6 Spin-charge separation 

One of the most famous and surprising results of one-dimensional many-body physics is the 
prediction of separate spin- and charge-excitations. We will show here that on the level of 
bosonization the spin-charge separation simply corresponds to taking suitable linear combina- 
tions of degenerate states. 

2.6.1 Spin and charge excitations 

All calculations that we have done so far actually also apply for electrons with spin if we use 
two species of bosons (p a with index cr =1,1, one for spin up and one for spin down. According 
to equation (|65p the bosons can immediately be related to fermionic densities with spin up 
and spin down. 

: fr R/L {x) : = -±=d x r R/L {x) a =|, I (104) 

If we are interested in the total charge density at a point x, we will use the sum of spin up 
and spin down 

: P C r/l( x ) ■ =■■ P R/L ( X ) : + : Pr/l( x ) ■ ( 105 ) 
Likewise, we can calculate the total spin density at point x by using the difference 

: Pr/l( x ) : =■ Pr/l( x ) : - : Pr/l( x ) : ( 106 ) 
Analogously, we can define new boson fields that correspond to the spin and charge densities 




(107) 



This defines new zero modes and oscillator mode operators (the index for left and right movers 
L/R is omitted here) 




(108) 



which have again canonical commutation relations analogous to (|5ip and (|64p as can easily be 
verified. The mode expansion of the fields (|107p in terms of the new operators (|108p is then 
the same as before. 

It is a straightforward exercise to write the bosonization formulas (|65p and (|79p in terms 
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of the new spin and charge fields, for example: 

ip R (x) oc exp ^i\f^n<t) R {x)^ = exp iV2ir ((f> R (x) + (p R (x)) 



: Pr( x ) : = -7^ d ^ R {x) = -7= (d x </>%(x) - d x <j) s R (x)) ( 109 ) 



: PR C ( X ) : = \[\ d x<t>R C { x ) 



Most importantly, the Hamiltonian (|66p can also be written in terms of the new spin and 
charge operators 

H = avp / dx ((d x <f> R ) 2 + (d x <p l R ) 2 + (d x( f>l) 2 + (d x cf> L ) 2 ) 

J ° " (110) 

= av F / dx ((d x <t> R ) 2 + {d x ^ s R ) 2 + {d^if + (d^i) 2 ) 
Jo 

Since the spin and charge excitations appear separately in the Hamiltonian (jllOp . the 
partition function factorizes. It is therefore possible to regard the spin and charge particles as 
the "new" independent fundamental excitations instead of the "old" spin up and spin down 
particles. Indeed, because of the degeneracy of the spin up and spin down channel, we could 
have used any canonical rotation to define new particles, but spin and charge are particularly 
useful when interactions are present. This is intuitively clear because realistic interactions 
will couple total charge densities and thereby lift the degeneracy between the spin and the 
charge channel. With interactions we are therefore forced to use the picture of spin and charge 
excitations, since the freedom of rotating degenerate channels is lost. 

It should be noted that the spin and charge separation is not obeyed exactly for the particle 
numbers, because the "new" zero modes still must obey the "old" quantization formula (|90p . 
In particular, the number of charge and spin particles is given by 

n C R = n R + n l R = J-Q C R 

V 1 (111) 

7T 



Both numbers must be integers, but cannot be changed independently. Adding a spin particle 
must always be accompanied by adding or removing a charge particle and vice versa, i.e. the 
total sum of spin and charge particles must always remain even. This just reflects the fact 
that we always have to add and remove real electrons, instead of spin/charge quasi-particles. 

We are now in the position to also interpret the nature of the spin and charge excitations. 
As discussed in section 2.2 we understand the bosonic excitations in terms of fermions being 
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t I t 4 




Figure 17: Charge excitations are symmetric linear combinations of up and down excitations. 
Spin excitations are anti-symmetric linear combinations of up and down. 

shifted up the spectrum. The spin and charge excitation now simply correspond to odd 
and even linear combinations of those shifts according to equation (|108j) . This is depicted 
schematically in diagram 1171 

2.6.2 Spin and charge correlation functions 

We now have the tools in order to calculate correlation functions of electrons with spin 

(^ Tt (x)V T (y)) « ae-^-f) (vI?(s)Vto)) +^ lkFix ~ v) (^Wi(y)) 

Using equation (|109p we can calculate the right moving part (using the same calculations as 
in section 2.5) 




(112) 



Without interactions the outcome is of course identical to the spinless result (197ft . but the 
calculation shows that the correlation function factorizes into a spin and a charge part. When 
interactions are introduced in the next section we expect that the degeneracy of spin and 
charge is lifted and therefore the correlations decay differently and independently for spin- 
and charge-like excitations. 
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3 Electron-electron interactions 



Finally we are now in the position to deal with interaction effects, which of course is the main 
goal of bosonization. After having derived the mathematical foundation of bosonization in 
the previous section, it is now straightforward to apply that prescription to typical interaction 
terms discussed above (see equations ([9]) and (fit))) ) 

3.1 Scattering processes 

Let us consider the standard model of density-density interactions in equation ([9]) 

N N 

Hint = E E ^( X i)^( X i) U ( m )^ (Xj+m)lp(x j+m ) 
j=l m=l 

or equivalently in terms of the wave-numbers k in equation (jlOp 

Hint = at E c\c k -.AkU{Ak)c\,c k > + Ak- (H3) 

k,k',Ak 

Following the standard program of bosonization, the first step is always to restrict the creation 
and annihilation operators in a linearized region \k-kp\<A and \k+kp\<A around kF according 
to Fig. O (see section 2.1). Let us systematically analyze the different possibilities as shown 
in Fig. ESI 

Forward scattering 

We first consider the case of small momentum transfer Ak < A, which is commonly referred 
to as forward scattering, since a right mover is always scattered into a right mover and a left 
mover is always scattered into a left mover. Using equation (|32p and restricting the momenta 
k and k' to left movers — kp — A < k, k' < —kp + A or right movers kF — A < k, k 1 < kp + A, 
the corresponding interaction Hamiltonian is given by 

-^forward = ^ {^k C k-Ak + C fc'^ c fc'-Afc) U {Ak) (cffc§ + & k + cJc k , +Ak ^j 

k,k',Ak 

1 X 

= Tr E {92 {p-AkPAk + P-AkPAk) + 94 (p-AkPAk + P-AkPAk)) 

iV Afc=-oo (114) 
1 00 

= A? E ( 2n ^ (&f Iff + «) + 2ng 4 (b*b* + btfbfy 

n=l 

+2g 2 n R n L + g 4 (n 2 R + n 2 L ) + const.) 
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Backward scattering Umklapp scattering 

Figure 18: The different scattering processes in the linearized approximation. 

where we have used equations f|41[) and (|50p . The constant comes from ordering the bosons 
and plays no significant role in the Hamiltonian. It is customary to denote the right-left 
interaction by gi =U(Ak), which is sometimes also called dispersion scattering. The density- 
density interaction strength on the same branch is denoted by =U(Ak). Often g 2 and 54 
are taken to be independent of momentum, which is justified for short range potentials U(m) 
that drop off quickly for m > 1/A. In that case, we can write the forward scattering term 
more compactly with the help of equation (|63p 

^forward = - dx (2g 2 (d x <f> R d x <t> L ) + <? 4 {(d x cf> R ) 2 + {d x <f> L ) 2 )) (115) 
71 Jo 

Backward scattering 

There are additional terms in equation (|10p corresponding to large momentum transfer 
around ±2kp, which scatter left into right movers and vice versa. We use a new parameter- 
ization AA; = ± (2kF + 5k), where 5k is now assumed to be small. Using (|lUp and (|32p . we 
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obtain 

^backward = c k^ c k~SkU {2k F + 6k)c$cg +Sk + c£_ Sk U (-2k F - 6k)t$ (%, +Sk 

k,k',Sk 

- 1 ( R] R L if , if I H «t\ 

— 31 \ C k C k+Ak c k'-Ak c k' + c fc C k+Ak c k'-Ak c k' J 

k,k',Ak 

where we have again re-parameterized Ak = 5k+k'—k in the last step and defined the backward 
scattering amplitude g\ = U (2k F + Ak + k — k'), which can be assumed to be approximately 
constant for k, k' and Ak on the range of the cutoff A in most cases. Therefore, using (|4ip 
and (1501) 

^backward ~ ^ 91 {PAkP-Ak + PAkP-Ak) (H6) 

k,k',Ak 

which has the same effect as the dispersion scattering equation ()114j) above and therefore can 
be absorbed right away by a redefinition of the corresponding scattering amplitude g2- For 
electrons with spin the backward scattering also introduces a spin-flip term, which can be 
taken into account by a renormalization group treatment, but will not be addressed here. 

Umklapp scattering 

Finally, there is the possibility of Umklapp scattering of two left movers into two right 
movers and vice versa in case that Ak = db {Akp + 5k) = ± (2ir + 5k). This Umklapp scattering 
(denoted by 53) is obviously only present if the special condition for half-filling 4k F = 2ir is 
fulfilled. The resulting operator can then also be bosonized with the help of the fermion field 
expressions (|79p and (|80p . but we will not consider this special case here. 

It should be noted here that the bosonization of the various scattering terms can also be 
derived from the interaction Hamiltonian in real space ([9]) by using the linearization formula 
(|34p together with (|65p if we assume a short range potential U{m) (as can be verified as 
an exercise). However, the derivation in momentum space has the advantage that the k- 
dependence of the scattering amplitudes 52 and c/4 can be preserved. 

3.2 The Luttinger Liquid parameter: Boguliobov transformation 

The bosonized interaction Hamiltonian in (|114p or (|1 15[) is much easier to treat than the 
fermionic version Q, because it is represented by terms that only involve a product of two 
boson operators (bilinear) as opposed to four fermion operators in (J9]) and (|10p . In fact, it is 
well known how to solve such problems by a suitable redefinition of the bosonic operators so 
that the interaction term has the same form as the original Hamiltonian, which only contains 
counting operators. This procedure is called a Boguliubov transformation and is given by the 
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general ansatz 

b 1 } = cosh 9 - b^smh.6 

(117) 

^ = b^ cosh 9 -65 sinhfl 



which defines new operators b% and 6£. It is left as an exercise to show that those new 
operators also obey the canonical commutation relations (|5ip . 

We notice that the 54 interaction in (|114p or (|115p is already in the non-interacting form 
of equations (|66p and (|67p and can be absorbed right away by defining a new Fermi velocity 

y F=Vp+ 9± (118) 
7T 

Therefore the complete Hamiltonian can be written as 

H = av F I dx ( {8 x (t> R ) 2 + {d x 4> L f + ^-d x </>Rd x <l> L ) 
Jo V v F ir ) 

= E ^ (*W + W + ^ ( W + «)) (H9) 



(where 52 is assumed to contain also any possible backward scattering contributions (|116p ). 
It is clear now that we want to use the canonical transformation (|117|) in order to get rid of 
the cross terms bn^bn^ + b^b^. This can be done by comparing with the following expression 
(up to a constant) 



bn ] bn + bffin = ( cosh2 + sinh2 e ) W^t + ^n h n) ~ 2 cosh sinh 0(6^6^ + b^) (120) 
which yields the following result for the rotation angle 9 



K = e^ = J 1 - 92/ ^ F7T (121) 
V 1 + 92/VF'n- 

where K is the so called Luttinger liquid parameter. It is straightforward to verify this result 
using the following formulas 

sinh.9 = -(y/K -l/VW^j cosh 9 = i( VK + l/VW^j sinh 9 + cosh 9 = \[K 

sinh 2 6» + cosh 2 (9 = (K + 1/K) /2 2 sinh 6» cosh 6» = (K - l/K) /2 (122) 
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Likewise, we can define new number operators 



Qr = Qr cosh 9 — Ql sinh 9 



(123) 



Ql = Ql cosh 9 — Qr sinh 9 



after which the Hamiltonian is in the standard form again 



CO 



2ttvf 



N 



n 



+ i L M) + v -§(qI + Ql) 



71=1 



(124) 



av F / dx \{d x (j) R ) 2 + (d x 4> L ) 2 



■Jo 



where we have also defined new fields 4>r and 4>L analogous to equations (I117P and (I123p . 
More compactly we can summarize these formulas in a canonical rescaling equation for the 
difference and the sum of the left and right moving fields 



This implies also that the creation operators 4>q and </>q are transformed, but it should be 
remembered that the quantization relation must still be fulfilled and particle excitations are 



In summary, we have therefore solved an interacting problem exactly by a simple canon- 
ical transformation. Moreover, all interactions can be described by a single Luttinger liquid 
parameter K = e 28 in equation (|12ip . Even though K can in principle be dependent on the 
momentum transfer k = n27r/N it is often sufficient to just use a constant value for short 
range interacting models. However, it turns out that the expression is only correct to low- 
est order in the actual scattering amplitudes g<i in typical lattice Hamiltonians like nearest 
neighbor interactions or the Hubbard model [16], because higher order operators renormalize 
the interactions. Therefore, the actual value of the Luttinger liquid parameter K must almost 
always be inferred from other theoretical methods [16,17] or by comparisons with experiments. 
For repulsive interactions we see from equation (|12ip that K < 1. For K = 1 we recover the 
non-interacting theory. 

3.3 Correlation functions 

To conclude this chapter, we would now like to apply our solution of the problem in order 
to calculate correlation functions in an interacting model. Again starting from (|9ip we are 





1 



(125) 



performed by the "old" operator exp (— z\/47T(/)^) . 
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interested in the correlator of the right moving field tp R oc j^*^ = e VMfecosh0+<^sinh<?) _ 
In particular, 

_ e i7rsign(a;-j/)/2 / e ~i^/In(<i> R (x)-<i> R (y)) cosh9\ / e ~i^/In(4> L (x)-4> L (y)) sinht 



where we have used (|79p and (|117p . The transformed operators <pn and 4>l have the same 
canonical mode expansion as before and therefore the calculation of the correlation functions 
proceeds as in section 2.5. 

In contrast to equation (|97p . however, the overall proportionality constant is now cut-off 
dependent and not normalized to unity. The left moving correlator differs by a minus sign. 



^1*0*0^(2/)) k i sign(x - y) 2£smj(x - y) 



- cosh 2 



TT 



21 sin — (x — y) 



- sinh 2 



oc i sign(a; — y) 2£sm—(x — y) 



TT 



-(K+l/K)/2 



(126) 



As in section 2.5 we can also use the time dependent mode expansion in order to calculate 
the Green's function. Accordingly, we find 



Vi(*)^fl(0)\ oc % sign(t) 



2£sin 



TTVpOt 



-(K+l/K)/2 



(127) 



A famous result for the density of states can now be obtained in the thermodynamic limit 
■ oo, since according to equations (fT7|) and (|127|) 



G R (t,x) = -i (Mx,t),tf(x,0)})0(t) t-^+W/ 2 



(128) 



and therefore from the density of state becomes, 



p(u) oc / dte^i-^+V^/a oc u {k+x/k)/2-i _ 



(129) 



where we have simply used a rescaling of the integration variable x= cut. This is the famous 
result of the depletion at low frequencies of the single particle spectral weight in a Luttinger 
liquid with a characteristic power-law 



p(uj) oc uj a , 



K + l/K-2 
a = - > 



(130) 



This depletion is one of the main hallmarks which is used in order to detect Luttinger liq- 
uid behavior experimentally with tunneling and photoemission experiments [7-9] as described 
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in section 1.2. It must be emphasized again that the proportionality constant is not known, 
but tunneling experiments cannot detect the overall amplitude of the density of states either. 
Theoretically it is possible to describe the proportionality constant of correlators phenomeno- 
logically with a cut-off energy scale [1-5] or by using a restricted momentum range for the 
interactions [18], but these are just mathematical tools that do not determine the actual value. 
In some microscopic models it is possible to fix the proportionality constant and the Luttinger 
liquid parameter by comparison with other exact methods [16,17,19]. 

As mentioned above, the correlators for electrons with spin factorize in a spin part and 
a charge part. In this case the interactions only transform the charge bosons unless spin- 
dependent scattering is considered. Therefore, we can immediately generalize equation (\112\) 
for the interacting case 

(V'fl (^LCy)) a i sign^ - y) 

and analogously for the time correlations. The corresponding result for the power-law depletion 
of the density of states becomes analogous to (|129p and (|130p 

p{uj) ocuj a , a = '- > (132) 

If a general correlation in space and time is calculated, the two factors in (I13ip contain different 
velocities v c > v s , corresponding to the separate excitations of spin and charge. 

This concludes the elementary introduction to bosonization. Using the tools we have 
developed here the reader is encouraged to explore also some of the more advanced topics, 
which can for example be found in the suggested review articles [1-5] . 



2^sin- 



(x + y) 



2£sin 



7T (X 



y) 



-(IX-t-L/FL )/1 



(131) 



4 Appendix: The boson cummulant formula 

The cummulant formula for bosons states that the finite temperature expectation value of 
an exponential of a linear combination of boson creation and annihilation operators can be 
expressed as the exponential of an expectation value in the following way 

^ e ab+/3b^ _ e ((ab+0f)/2 ( 133 ) 

If more than one species of bosons is present, this expression factorizes on both sides so that it 
can be used for any linear combination of bosons f = ^2 [ct n bn + flnbn) • It is not valid for the 
zero-mode operators, however (see Ref. [15] for a discussion on finite temperature expectation 
values of zero modes). 
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The right hand side of equation f)133|) can be evaluated rather easily by use of the Bose- 
Einstein distribution for finite temperatures /3=1//cbT in terms of the corresponding energy 
quantum Hu 

oo 

j2 e- 71 ^ (n\tfb\n) 



g— n/3fru) 

n=0 



e -f3Hui _ l 



Therefore, 



e ((ab+0) 2 )/2 _ e a/3(66t+6t6)/2 

= e 1 -" 2 



(134) 



where we have introduced the Boltzmann weight q = exp(-/3hu>). 

For the left hand side of (|133p we use the Baker Hausdorf formula (|82p and expand the 
exponential 



-a/3/2 



E 

n,n'=0 



a n /3 n 
n\n'\ 



-a/3/2 \ ^ m=0 

n=0 1 j £ g m 

m=0 

n=0 m=0 



m ; 



(135) 



6 n 6 n t 



77) 



where we have used the standard expression for a temperature expectation value with the 
Boltzmann weight q. For a harmonic oscillator it is well known that the expectation value in 
the m th excited state is given by 
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h n h n\ 



m 



ml 



(e.g. by repeated use of $ \m) = \Jm + 1 \m + 1) ). Finally, we insert into equation (|135H the 
following Taylor expansion around q = 



a-^-^E 



m=0 



(m + n)\ 
n\m\ 
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(which can be verified by repeated differentiation). Therefore, equation fj 135|) becomes 



,aH-/96t\ = (1 _ q)e - a /3 /2 - q)~n-t 

n=0 n ' 

Z-, n! ^1-^ 



e 

n=0 

a/3 



e a ^^ 2 exp 



1 



which is exactly the same expression as on the left hand side (|134p and therefore concludes 
the proof of the cummulant theorem. Note that at T = 0, we have q = exp(-/?foi;)— > 0, so 
that 

= e ((a6+/36t)2)/ 2 = gxp 

as already shown in equation ([5 
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